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Abstract 


This document contains graphical data necessary for the preliminar\ design of ballistic 
missions to Venus Contours of launch energy requirements, as \^ell as many other launch 
and Venus arrival parameters, are presented in launch date arm aJ date space for all 
launch opportunities from 1991 through 2005. In addition, an e.xtensive text is included 
which explains mission design methods, from launch window development to Venus 
probe and orbiter arrival design, utilizing the graphical data in this volume as well as 
numerous equations relating various parameters. This is one of a planned series of mission 
design documents which will applv to ai! planets and some other bodies in the solar 
b\ stem 
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Preface 

This publication is one of a series of volumes devoted to interplanetary trajectories of 
different types. Volume 1 deals with ballistic trajectories. The present publication is 
Part 1 and describes ballistic trajectories to Venus, Part 3, which was published in 1982, 
treated ballistic trajectories to Jupiter. Parts 2 and 4, which were published earlier in 
1983, described ballistic trajectories to Mars and Saturn, respectively. Part 5, which will 
be published in the near future, will treat ballistic Mars-to-Earth return trajectoiies. 
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I. Introduction 

The purpose of this series of Mission Design Handbooks is 
to provide trajectory designers and mission planners with 
graphical trajectory information, sufficient for preliminary 
interplanetary mission design and evaluation. In most respects 
the series is a continuation of the previous three volumes of 
the Mission Design Data, TM 33-736 (Ref. 1 ) and its predeces- 
sors (e.g.. Ref. 2); it extends their coverage to departures 
through the year 2005 A.D. 

The entire series is planned as a sequence of volumes, each 
describing a distinct mission mode as follows: 

Volume I: Ballistic (i.e., unpowered) transfers be- 

tween Earth and a planet, consisting of 
one-leg trajectory arcs. For Venus and 
Man missions the planet-to-Earth return 
trajectory data are also provided. 

Volume II: Gravity-Assist (G/A) trajectory transfers, 
comprising from two to four ballistic 
interplanetary legs, connected by suc- 
cessive planetary swingbys. 

Volume III: Delta- V-EGA (AVEGA) transfer trajec- 
tories utilizing an impulsive deep-space 
phasing and shaping bum, followed by a 
return to Earth for a G/A swingby maneu- 
ver taking the :>pacecraft (S/C) to the 
eventual target planet. 

Each volume consists of several parts, describing trajectory 
opportunities for missions toward specific target or swingby 
bodies. 

This Volume I. Part 1 of the series is devoted to ballistic 
transfen between Earth and Venus. It describes trajectories 
taking from 100 to 500 days of flight time for the 10 succes- 
sive mission opportunities, departing Earth in the following 
years. 1991, 1992/3. 1994, 1996. 1997/8. 1999, 2000/1, 
2002, 2004, and 2005/6. 

Individual variables presented herein are described in detail 
in subsequent sections and summarized again in Section IV. 
Suffice it to say here that all the data are presented in sets of 
11 contour plots each, displayed on the launch date/arrival 
date space for each opportunity. Required departure energy 
C 3 , departure asymptote declination and right ascension, 
arrival and its equatorial directions, as well as Sun and 
Earth direction angles with respect to the departure/arrival 
asymptotes, are presented. 

It should be noted that parts of the launch space covered 
may require launcher energies not presently (1983) available, 


but certainly not unrealistic using future orbital assembly 
techniques. 

A separate series of volumes (Ref. 3) is being published 
concunently to provide purely geometrical (i.e., trajectory- 
independent) data on planetary positions and viewing/orienta- 
tion angles, experienced by a spacecraft in the vicinity of these 
planetary bodies. The data cover the time span through 2020 
A.D.. in order to allow sufficient mission duration time for all 
Earth departures, up to 2005 A.D. 

Tlte geometric data are presented in graphical form and 
consist of 26 quantities, combined into eight plots for each 
calendar year and each target planet. The graphs display equa- 
torial declination and right ascension of Earth an^ Sun (plan- 
etocentric), as well as those of the target planet (geocent ic); 
heliocentric (ecliptic) longitude of the planet, its heliocentric 
and geocentric distance; cone angles of Earth and Canopus, 
clock angle of Earth (when Sun/Canopus-oriented); Earth-Sun- 
planet, as well as Sun-Earth-planet angles; and finaUy, rise and 
set times for six deep-space tracking stations assuming a 6 -deg 
horizon mask. This information is similar to that in the second 
part of each of the volumes previously published (Ref. 1 ). 

II. Computational Algorithms 

A. General Description 

The plots for the entire series were computer-generated. A 
minimum of editorial and graphic support was postulated from 
the outset in an effort to reduce cost. 

A number of computer programs were created and/or modi- 
fied to suit the needs of the Handbook production. 

The computing effort involved the generation of arrays of 
transfer trajectory arcs connecting departure and arrival 
planets on a large number of suitable dates at each ^ody. 
Algorithms (computational models) to solve this problem can 
vary greatly as to their complexity, cost of data generated, and 
resulting data accuracy. In light of these considerations, the 
choice of n^ethods used in this effort has been assessed. 

B. Two-Body Conic Transfer 

Each departure/arrival date combination represents a unique 
transfer trajectory between two specified bodies, if the number 
of revolutions of the spacecraft about the primary (e.g,, the 
Sun) is specified. The Lambert Theorem provides a suitable 
framework for the computation of such primary-centered 
trajectories, but it is of practical usefulness only if restricted 
two-body conic motion prevails. 

Restricted two-body motion implies that the dynamical 
system consists of only two bodies, one of which, the primary. 
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is so much more massive than the other, that all of the system’s 
gravitational attraction may be assumed as concentrated at a 
point-the center of that primary body. The secondary body 
of negligible mass (e.g., the spaceciaft) then moves in Kep- 
lerian (conic) orbits about the primary (e.g., the Sun) in such 
a way that the center of the primary is located at one of the 
foci of the conic (an ellipse, parabola, or hyperbola). 

The Lambert Theorem states that given a value of the gravi- 
tational parameter (j. (also known as GM) for the central body, 
the time of flight between two arbitrary points in space, 
Rj and R 2 * is a function of only three independent variables: 
the sum of the distances of the two points from the focus, 
IRj I + IR 2 I, the distance between the two points C= |Rj - Rj I, 
and the semimajor axis, a, of the conic orbital flight path 
between them (Fig. 1 .) 



Detailed algorithm descriptions of the Lambert method, in- 
cluding necessary branching and singularity precautions, are 
presented in numerous publications, e.g.. Refs. 2 and 4. The 
computations result in a set of conic classical elements {a, e, /. 
n, cu, Uj) and the transfer angle, AUj 2 » or two equivalent 
spacecraft heliocentric velocity vectors, depar- 

ture, the other at the arrival planet. Subtraction of the appro- 
priate planetary heliocentric velocity vector, at 

the two corresponding times from each of these two space- 
craft velocity vectors results in a pair of planetocentric velocity 
states “at infinity” with respect to each planet (Fig. 2): 

I ^5/Cati ^PLANET i ^ 

where 1 = 1 and 2 refer to positions at departure and arrival, 
respectively. The scalar of th^s vector is also referred to as 


Fig. 1. The Lambert problem geometry 


the hyperbolic excess velocity, “K-infinity” or simply “speed” 
(e.g.. Ref. 4). The represents the velocity of the spacecraft 
at a great distance from the planet (where its gravitational 
attraction is practically negligible). It is attained when the 
spacecraft has climbed away from the departure planet, fol- 
lowing injection at velocity V^: 



or before it starts its fall into the arrival planet’s gravity well, 
where it eventually reaches a closest approach (periapse, C/4) 
velocity 


Vi, 

^FARTH 



Fig. 2. Departure geometry and velocity vector diagram 
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The variables Tj and refer to the departure injection and 
arrival periapse planetocentric radii, respectively. Values for 
the gravitational parameter ix (or GM) are given in subsequent 
Section V on constants. 

The vectors, computed by the Lambert method, repre- 
sent a body center to body center transfer. They can, however, 
be translated parallel to themselves at either body without 
excessive error due to the offset, and a great variety of realistic 
departure and arrival trajectories may thus be constructed 
through their use, to be discussed later. The magnitude and 
direction of as well as the angles that this vector forms 
with the Sun and Earth direction vectors at each terminus are 
required for these mission design exercises. 

Missions to the r *latively small terrestrial planets such as 
Venus are suited to be analyzed by the Lambert nrnhod. as 
the problem can be adequately represented by the restricted 
two-body formulation, resulting in flight time errors of less 
than 1 day - an accuracy that cannot even be read from the 
contour plots presented in this document. 

C. Pseudostate Method 

Actual precision interplanetary transfer trajectories, espe- 
cially those involving the giant outer planets, do noticeably 
violate the assumptions inherent in the Lambert Theorem. 
The restricted two-body problem, on which that theorem is 
based, is supposed to describe the conic motion of a massless 
secondary (i.e., the spacecraft) about the point mass of a pri- 
mary attractive body (i.e., the Sun), both objects being placed 
in an otherwise empty Universe. In reality, the gravitational 
attraction of either departure or target body may significantly 
alter the entire transfer trajectory. 

Numerical A-body trajectory integration could be called 
upon to represent the true physical model for the laws of 
motion, but would be too costly, considering the number of 
complete trajectories required to fully search and describe a 
given mission opportunity. 

The pseudostate theory, first introduced by S. W. Wilson 
(Ref. 5) and modified to solve the three-body Lambert prob- 
lem by D. V. Byrnes (Ref. 6), represents an extremely useful 
improvement over the standard Lambert solution. For the 
giant planet missions, it can correct about 95 percent of the 
three-body errors incurred, e.g., up to 30 days in flight time on 
a typical Jupiter-bound journey. 


Pseudostate theory is based on the assumption that for 
modest gravitational perturbations the spacecraft conic motion 
about the primary anu the pseudo-conic displacement due to a 
third body may be superimposed, if certain rules are followed. 

The method, as applied to transfer trajectory generati^' 
does not provide a flight path -only its end states. It solv 
original Lambert problem, however, not De.ween the le 
planetary positions themselves, but instead, between two • im- 
puted "pseudostates." These are obtained by iteration on twe 
displacement vectors off the planetary ephemeris positions on 
the dates of departure and arrival. By a suitable superposition 
with a planetocentric rectilinear unnact hyperbola and a 
constant-velocity, "zero gravity," sweepback at each end of 
the Lambertian conic (see Fig. 3), a satisfactory match is 
obtained. 

Of the five arcs involved in the iteration, the last three 
(towards and at the target planet) act over the full flight time. 
ATjj. and represent: 

(1) Conic heliocentric motion between the two pseudo- 
states R* and R* (capital R is used here for all helio- 
centric positions), 

(2) The transformation of R* to a planetocentric position, 

(lower case r is used for planetocentric positions), 
performed in the usual manner is followed by a 
“constant -velocity” sweepback in time to a point 
r* = rj - V^2 ^ ^^12’ correcting the planetocentric 
position to what it would have been at , had 
there been no solar attraction during and 

finally, 

(3) The planetocentric rectilinear incoming hyperbola, 
characterized by: incoming F-infuiity V^^, a radial 


Fig. 3. Fseudostate transfer geometry 
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target planet impact, and a trip-time ATjj from r* 
to periapse, which can be satisfied by iteration on the 
r^’inagnitude and thus also on R^. 

This last aspect provides for a great simplification of the 
formulation as the R* end-point locus now moves only along 
the Vooj II Fj vector direction. The resulting reduction in com- 
puting cost is significant, and the equivalence to the Lambert 
point-to-point conic transfer model is attractive. 

The first two segments of the transfer associated with the 
departure planet may be treated in a like marmer. If the planet 
is Earth, the pseudostate correction may oe disregarded (i.e., 
R, =R,), or else the duration of Earth's perturbative effect 
may be reduced to a fraction of 

equal a fixed quantity, e.g., AT. = 20 days. The latter value 
was in fact used at the Earth’s side of the transfer in the 
data generation process for this document. 

The rectilinear pseudustate method described above and in 
Ref. 7 thus involves an iterative procedure, utilizing the 
standard Lambert algorithm to obtain a starting set of values 
for at each end of the transfer arc. This first guess is then 
improved by allowing the planetocenlric pseudostate position 
vector to be scaled up and down, using a suitable partial at 
either body, such that the time required to fall along 

the rectilinear hyperbola through (the sum of th<^ sweep- 
back distance X AT} and the planetocenlric distdiice jr,* |), 
equal the gravitational perturbation duration, AT}. Both and 
along which the rectilinear fall occurs, are continuously 
leset utilizing the latest values of magnitude and direction of 
at each end, /, of the new Lambert transfer arc, as the 
iteration progresses. The procedure converges rapidly as the 
hyperbolic trip time discrepancy, AAT= ATj^^q - AT}, falls 
below a preset small tolerance. Once the vectors at each 
planet are converged upon, the desired output variables can be 
generated and contour plotted by existing standard algorithms. 

As mentioned before, the pseudor.tate method was found 
to be unnecessary for the accuracy of this Venus-oriented 
handbook. The simpler Lambert method was used instead in 
the dala-generation process. 

III. Trajectory Characteristics' 

A. Mission Space 

All realistic launch and injection vehicles are energy -limited 
and impose very stringent constraints on the interplanetary 
mission selection process. Only those transfer opportunities 
which occur near the times of a minimum Earth departure 
energy requirement are thus of practical interest. On either 
side of such an optimal date, departure energy UiCre^ses, first 
slowly, followed by a rapid increase, thus requirmg either a 
greater launch capability, or altcrnativelv a lower allowable 


payload mass. A **launch period,” measured in days or even 
weeks, is thus definable: on any day within its confines the 
capability of a given launch/injecticn vehicle must equal or 
exceeu the departure energy requirement for a specified 
payload weight. 

In the course of time these minimum departure energy 
opportunities do recur regularly, at “synodic period” intervals, 
reflecting a repetition of the relative angular geometry of the 
two planets. If cu, and are the orbital angular ra.es of the 
inner and outer of the two planets, respectively, moving about 
the Sun in circular orbits, then the mutual configuration of 
the two bodies changes at the following rat^ 

^12 ’ " ^2 ’ 

If a period of revolution, Py is defined as 


then 


where the synodic period, is the period of planetary 
geometry recurrence, while and P 2 are the orbital “sidereal 
(i. .., inertial) periods” of the inner (faster) and the outer 
(slower) planet considered, respectively. 

Since planetary orbits are neither exactly circular nor cr 
planar, launch opportunities do not repeat exactly, 
years being better than others in energy requirements o. . 
other parameters. A complete repeat of trajectory chameter 
istics occurs only when exactly the same orbital geometry of 
departure and arrival body recurs. For negligibly perturbed 
planets approximately identical inertial positions in space at 
departure and arrival imply near-recurrence of transfer tra- 
jectory characteristics. Such events can rigorously be assessed 
only for nearly resonant nonprecessing planetary orbi> , it., 
for those whose periods can be related in terms of uteg<*r 
fractions. For instance, if five revolutions of one body c. ' 
respond to three revolutions of the other, that time interval 
would constitute the “period of repeated characteristics.” 
Near-integer ratios provide nearly repetitive configurations 
with respect to the lines of apsides and nodes. The Earth- 
relative synodic period of Venus is 583.920 days, i.e., about 
1.599 years. Each cycle of 5 consecutive Venusian mission 
opnortunities amounts to 2919.602 days and is nearly repeti- 
tive. driven by 13 Venusian sidereal periods of 224.7007 days. 
It is obvious that for an identical mutual angular geometry 
Venus would be found short of its inertial position in the pre- 
vious cycle by 1.51 days worth of motion (2.414 deg), while 
Earth would have completed 2919.602/365.25 = 7.9934 revo- 
lutions, being short of the old mark by the same angular 
amount as Venus. 
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Repeated occurrences of this e icellent S-year. 5-oppor- 
tunity cycl*^ do not provide for an improved geometric con- 
figuration: in fact, a slow deterioration of repeatability is 
observed. 

A variety of considerations force the realistic launch period 
not to occur at the minimum energy combination of depar- 
ture and arrival dates. Launch vehicle readiness status, proce- 
dure slippage, weather anomalies, multiple launch strategies, 
arrival characteristics-all cause ihe launch or, more generally, 
the departure period to be extended over a number of days or 
weeks and not necessarily centered on the minimum energy 
date. 

For tills document, a 200-day departure dat^ coverage span 
was selected, primarily in order ewcompass launch energy 


requirements of up to a C 3 = 30 kni^/s^ contour, where 
C 3 = i.e.. twice the injection energy per unit mass. 

Ej = Arrival date coverage was set at 180 days to 

display *n Mons from 75-250 days of flight time. 

The matrix of depai.vre and arrival dates to be presented 
comprises the “mission space” for each departure opportunity. 

B. Transfer Tra|ectory 

As previously stated, each pair of departure/arrival dates 
specific, unique transfer trajectory. Each such point in the 
mission space has associated with it an array of descriptive 
variables. Departure energy, characterized by C3, is by far the 
most significant among these parameters. It increases towards 
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the edges of the mission space, but it also experiences a dra- 
matic rise along a “ridge,” passing diagonally from lower left 
to upper right across the mission space (Fig. 4). This distur- 
bance is associated with all diametric, i.e., near-180-deg, 
transfer trajectories (Fig. 5). 

In 3-dimensionaJ space the fact that all planetary orbits are 
not strictly coplanar causes such diametric transfer arcs to 
require high ecliptic inclinations, culminating in a polar flight 
path for an exact 180-deg ecliptic longitude increment be- 
tween departure and arrival points. The reason for this 
behavior is, as shown in Fig. 5, tliat the Sun and both trajec- 
tory end points must lie in a single plane, while they are also 
lining up along the same diameter across the ecliptic. The 
slightest target planet orbital inclination causes a deviation 

V = TRUE TRANSFER 
ANGLE 

- TRANSFER ANGLE 
POLAR IN ECLIPTIC PLANE 



Fig. 5. Effect of transfer angle upon inclination of 
trajectory arc 



Fig. 6. Nodal transfer geometry 


out of the ecliptic and forces a polar 180-deg transfer, in order 
to pick up the target’s vertical out-of-plane displacement. 

The obvious scle exception to this rule is the nodal transfer 
mission, where departure occurs at one node of the carget 
planet orbit plane with the ecliptic, whereas arrival occurs at 
the opposite such node. In these special cases, which recur 
every half of the repeatability cycle, discussed in the preceding 
paragraph, the transfer trajectory plane is indeterminate and 
may as well lie in the departure planet’s orbit plane, thus 
requiring a lesser departure energy (Fig. 6). The opposite strat- 
egy (i.e.. a transfer m (he arrival planet's orbit plane) may be 
preferred if arrival enerey, Voo,, is to be minimized (Fig, 7). 

It should be noted that nodal transfers, being associaied 
with a specitlc Venus arrival date, may show up in the data on 
several consecutive Venus miss.on opportunity graphs (e.g.. 
1996-2004). al corresponding to the paiticular nodal 

arrivi;! date. Their mission space position moves, from oppor- 
tunity to opportunity, along the 180-deg transfer ridge, b> 
gradually sliding towards shorter trip times and earlier relative 
departure dates. Only one of these opportunities would occur 
at or near the minimum departure ent.gy or the minimum 
arrival date, which require a near-pe ihelion to near- 
aphelion transfer trajectory. These pseudoT!ohmann nodal 
transfer opportunities provide significant energy advantages, 
but represent singularities, i.e., single-time-pomt missions, 
with extremely high error sensitivities. Present-dav mission 
planning does not allow single fixed-time departure strategies: 
however, future operations modes, eg., space station “on- 
time" launch, or alternately Earth gravity assist ( repeat. d) 
encounter at a specific time, may allow the advantages of a 
nodal transfer to be utilized in full. 

The T deg transfer ridge subdivides the mission space 
into two basic regions: the Type I trajectory space below the 
ridge, exhibiting less than 180-deg transfer arcs, and the 
Type II space whose transfers are longer than 180 deg. In 
general the first type also provides shorter trip times. 

Trajectories of both types are further subdivided in two 
parts-Classes 1 and 2. These are separated, generally hori- 
zontallv, by a boundary representing the locus of lowest 
energy for each departure date. Classes separate longer 
duration missions from shorter ones within each type. Type I, 
Class I missions could thus be prefer' :d becau:;e of their 
shorter trip times. 

Transfer energies become extremely high for very short 
trip times, infinite if launch date equals arrival date, and of 
course, meaningless for negative trip times. 

The reason that high-inclination transfers, as found along 
the ridge, also require such high energy expenditures at depar- 
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Fig. 7. Mission space with nodal transfer 


ture * a the spacecraft velocity vector due to the Earth's 
orbital velocity must be rotateu through large angles out of the 
ecliptic in addition to the need to acquire the required transfer 
trajectory energy. The value of C3 on the ridge is large but 
finite; its saddle point minimum value occurs for a pseudo- 
Hohmann fi.e., perihelion to aphelion) polar transfer, requiring 

where - 29.766 km;s. the Earth’s heliocentric orbital 
velocity, and = .723 AU. Venus* (the arrival planet's) semi- 
major axis. By a similar estimate, it can be shown that for a 
true nodal pseudo-Hohi.iann transfer, the minimum energy 
required would reduce to 


2 

=5^ 6 km^/s^ 

( 8 ) 

This is the lowest value of C 3 required to fly from Earth to 
Venus, assuming circular coplanar planetary orbits. 

Arrival K-infinity, . is at its lowest when the transfer 
trajectory is near-coplanar and tangential to the target planet 
orbit at arrival. 

Both C3 and Voo^ near the ridge can be significantly low- 
ered if deep-space deterministic maneuvers are introduced into 
the mission. The *‘broken-plane" maneuvers are a category of 
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such ridge-counteractiiig measures, which can nearly eliminate 
all vestiges of the near- 1 80-deg transfer difficulties. 

The basic principle employed in broken-plane transfers is to 
avoid high ecliptic inclinations of the trajectory by performing 
a plane change maneuver in the general vicinity of the halfway 
point, y:d\ that it would correct the spacecraft's aim toward 
the target planet’s oUt-of-ecliptic position (Fig. 8 ). 

Graphical data can be presented for this type of mission, 
but it requires an optimization of the sum of critical AT 
expenditures. The decision on which ^Ks should be included 
must be based on some knowledge of overall staging and arrival 
intentions, e.g., departure injection and arrival orbit insertion 
vehicle capabilities and geometric constraints or objectives 
contemplated. As an illustration, a sketch of resulting contours 
of C 3 is shown in Fig. 9 for a lypical broken-plane oppor- 
tunity represented as a narrow strip covering the ridge area 
(Class 2 of Type 1 and Class 1 of Type II) on a nominal 1991 


POLAR 

TRANSFER 



Fig. 8. Broken-plane transfer geometry 
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Fig. 9. Sketch of mission space with broken-plane transfer, effective energy requirements 
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launch/arrival date C3 contour plot. The deep-space maneuver 
is transformed into a C3 equivalent by converting the 
new broken-plane to an injection velocity at parking 

orbit altitude, adding and converting the sum to a new 

aiid slightly larger C3 value at each point on the strip. 

C. Launch/lnjection Geometry 

The primary problem in departure trajectory design is to 
match the mission-required outgoing K-infinity vector, V^, to 
the specified launch site location on the rotating Earth. The 
site is defined by its geocentric latitude, 0^, and geographic 
east longitude, (Figs. 10 and 1 1). 

Range safety considerations prohibit overflight of popu- 
lated or coastal areas by the ascending launch vehicle. For each 
launch site (e.g., Kennedy Space Center, Western Test Range, 
or Guiana Space Center), a sector of allowed azimuth firing 
directions is defined (measured in the site’s local horizontal 
plane, clockwise from north). For each launch vehicle, the 
allowed sector may be further constrained by other safety 
considerations, such as spent stage impact locations down the 
range and/or aown-range significant event tracking capabilities. 

The outgoing r-infinity vector is a slowly vary ing function 
of departure and arrival date and may be considered constant 
for a given day of launch. It is usually specified by its energy 


magnitude C3 = IV^p, called out as C^L in the plots and 
representing twice the kinetic energy (per kilogram of injected 
mass) which must be matched by launch vc cle capabilities, 
and the T-infinity direction with respect to the inertial Earth 
Mean Equator and equinox of 1950.0 (EME5U) coordinate 
system: the declination (i.e., latitude) of the outgoing asymp- 
tote 5^ (called DLA), and its right ascension (i.e., equatorial 
east longitude from vernal equinox, T) (or RLA). These 
three quantities are contour-plotted in the handbook data pre- 
sented in this volume. 

1. Launch azimuth problem. The first requirement to be 
met by the trajectory analyst is to establish the orientation of 
the ascent trajectory plane (Ref. 8). In its simplest form this 
plane rm»st contain the outgoing V-infinity (DLA, RLAl 
vector, the center of Earth, titui the launch site at lift-off 
(Fig. 10). As the launch site partakes ii. the sidereal rotation of 
the Earth, the continuously changing ascent plane manifests 
itself m a monotonic increase of the launch azimuth, with 
lift-off lime, (or its angUi,.: counterpart, measured 
in the equator plane): 

cos 0^ X tan 5^ - sin 0^ X cos (a^ - Q^) 


(9) 


ASYMPTOTE 



Fig. 10. Launch/injection trajectory plane geometry 
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The daJy time history of azimuth can be obtained from Eq. 
(9) for 1 given departure asymptote direction by 

following quadrant rules explained below and using the 
following approximate expressions (see Fig. II): 


^ ^^^DATE ^ EARTH ^ 


a 1 


DATE 


= 100.075 + 0.9856123008 X d 


so 


( 10 ) 

( 11 ) 


w ;ere 


Gii I 


DATE 


^EARTH 


'SO 


Right ascension of launch site (0^* at 
(deg) 

Greenwich hour angle at 0^ GMT of any 
date, the eastward angle between vernal 
equinox and the Greenwich meridian (deg), 
assumes equator is EME50.0 

sidereal (inertial) rotation rate of Earth 
(15.041067179 deg/Ti of mean solar time) 

launch date in terms of full integer days 
elapsed since 0^ Jan. 1 . 1950 (days) 


GREENWICH 
MERIDIAN 
0^^ GMT 

EAST LONGITUDE 
LAUNCH SITE \ 


LAUNCH SITE 
AT 0^ GMT 


ANGLE TRAVERSED 
BY LAUNCH SITE 
SINCE 0^ GMT 
OF LAUNCH DATE 


GREENWICH HOUR ANGLE 
(GHA) 





VERNAL 

EQUINOX 


DEPARTURE 

ASYMPTOTE 

MERIDIAN 


LAUNCH SITE 
AT LAUNCH, fGMT) 


Fig. 11. Earth equator plane definition oi angles involved 
in the launch problem 


tive) and 360.0 deg in the fourth quadrant (when cotan 2 
is negative). 


= lift-off time (h, GMT. i.e., mean solar lime) A generalized plot of relative launch time vs launch 

azimuth 2, can be constructed based on Eqs. (9) and (12), 


A relative launch time, measured with respect to an 

in rtial reference (the departure asymptote meridian's right 
ascention c<^). can be defined as a sidereal time (Earth’s 
ro‘ uion rate is 15.0 deg/h of sidereal time, exactly): 


w = T5:r 


( 12 ) 


This time represents a generalized sidereal time of launch, 
elapsed since the launch site last passed the departure asymp- 
tote niendian. a . 


if a fixed launch site latitude is adopted (e.g., 0j^ = 28.3 deg 
for Kennedy Space Flight Center at Cape Canaveral, Florida). 
Such a plot is presented in Fig. 12 with departure asymptote 
declination 6^ as the contour parameter. The plot is applicable 
to any realistic departure condition, independent of , date, 
or true launch time (Ref. 9). 

2. Daily launch windows. Inspection of Fig. 12 indicates 
that generally two contours exist for each declination value 
(e g-, = -10 deg), one occurring at during the a.m. 

hours, the other in the p.m. hours of the asymptote relative 
“day.” 


The cCtuil Greenwich Mea.i (solar) Time (GMT) of launch, 
may be obtained fi 'rr* r^^y by transforming it to mean 
solar lime and adding a uate. site, and asymptote-dependent 
adjustment: 

',«XI5.0 a (EAST) 

■■ ^ 

^EARTH ^EARTH 

'T expre'uon for 2^ (Eq. 9) must be used with computa- 
tional regaid for quadranU, singular points, and sign conven- 
tions. If the launch is known to be direct (i.e., eastward), 
then 2^ , when c ’4an 2^ is negative, must be corrected to 
2^ = 2/ ^ 180.0. For retrograde (westward) launches, 

180.0 de;^ must be added in the third (when cotan 2 is posi- 


Since lift-off times are bounded by preselected launch-site- 
dependent limiting values of launch azimuth 2^ (e.g.. 70 deg 
and 1 1 5 deg), each of the two declination contours thus con- 
tains a segment during which launch is permissible -“a launch 
window.” The two segments on the plot do define the two 
available daily launch windows. 

As can be seen from Fig. 12, for = 0, the two daily 
launch opportunities arc separated by exactly 12 hours; with 
an increasing I they close in on each other, until at 1 ~ 
|0y| they merge into a single daily opportunity. For |6^| > 
10y|, a “split” of that single launch window occurs, disallow- 
ing an ever-increasing sector of azimuth values. This sector is 
symmetric about east and its limits can be determined from 
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Fig. 12. Generalized relative launch time v$ launch azimuth and departure asymptote decimation 5^. 
Pair of typical example launch windows for 5^ = -10 deg shown by bold curve segments 

(reproduced from Ref. 9). 


sin 


V 


COS 6^ 

cos 


04) 


As 6^ gets longer, the sector of unavailable launch azimuths 
reaches the safety boundaries of permissible launches, and 
planar launch ceases to exist (Fig. 13). This subject will be 
addressed again in the discussion of “dogleg" ascents. 


Figure 14 is a sketch of a typical daily launch geometry 
situation, sh 'wn upon a Mercator map of the celestial sphere. 
The two launch windows exhibit a similar geometry since the 
inclinations of the ascent trajectory planes are functions of 
launch site latitude <t>i and azimuth Zj only: 


cos / = cos <t>j X sin 


(15) 


The angular equatorial distance between the ascending node 
and the launch site meridian is given by 

sin 0, X sin 2, 

sin (q, - n) = ^^- 7 —— (16) 

Quadrant rules for this equation involve the observation that a 
negative cos 2^ places (ct^ - H) into the second or third quad- 
rant. while the sign of sin (a, - H) determines the choice 
between them. 

The range angle 0 is measured in the inertial ascent trajec- 
tory plane from the lift-off point at launch all the way to the 
departure asymptote direction, and can be computed for a 
given launch time (j or azimuth 2^ 

already known from Eq. (9) as follows: 


cos d = sin 6^ X sin 0^ + cos 6^ X cos 0^ X cos (o^ ~ cHj) 


(17) 


The two daily opportunities do differ greatly, however, in 
the right ascension of the ascending mode of the orbit and 
in the length of the traversed in-plane arc, the range angle 0. 
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-90 -70 -50 -30 -10 10 30 50 70 90 

ASYMPTOTE DECLINATION . deg 


Fig. 13. Permissible regions of azimuth vs asymptote 
declination launch space for Cape Canaveral 


Sind = 


■ ®/.)X cos6„ 


(18) 


The extent of range angle 6 can be anywhere between 
0 deg and 360 deg, so both cos $ and sin $ may be desired in 
its determination. The range angle 6 is related to the equa- 
torial plane angle, da = ^ <l*scussed before Even 

though the two angles are measured in different planes, they 
both represent the angular distance between launch and depar- 
ture asymptote, and hence they traverse the same number of 
quadrants. 


Figure 1 5 represents a generaUy apphcable plot of centra] 
range angle 6 vs the departure asymptote declination and 
launch azimuth, computed using Eqs. (17) and (18) and a 
launch site latitude 0^ = 28.3 (Cape Canaveral). The twin 
daily launch opportunities are again evident, showing the 
significant difference in avaUable range angle when following a 
vertical, constant 6«, line. 


It is sometimes convenient to reverse the computational 
procedure and determine launch azimuth from known range 
angle $ and i.e., da = as follows: 


ASYMPTOTE RELATIVE LAUNCH TIME 
3 6 9 12 15 


O) 

m 

T3 

LU 

O 

D 


QC 

O 

H 

< 

D 

O 

LU 



90 135 180 225 

RIGHT ASCENSION a, deg 
(EARTH EQUATORIAL CELZSTIAL LONGITUDE) 


270 


315 


BY DEFINITION 

0^ = 24^ AT 
^-MERIDIAN 


DAILY TRACK 
OF LAUNCH SITE 
AT 0^ = 28.3® 


ASSUMED 
ASYMPTOTE 
LOCATION: 
6 <«=- 15 ® 
t^ = 315® 


Fig. 14. Typical launch geometry example in celestial (inertial) Mercator coordinates 
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400 



-70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 


DEPARTURE ASYMPTOTE DECLINATION . deg 

Fig. 15. Central range angle 0 between launch site and outgoing asymptote direction vs its declination and launch azimuth. 
Pair of typical example launch windows for = -10 deg shown by bold line segments 

(reproduced from Ref. 8). 


sin 2^ 


cos 6^ X sin (a^ - a^) 
sin 6 


sin 6^ “ cos 0 X sin 0^ 

cos X, = : — 7 

^ sin 0 X cos 0^ 


with the lighting conditions at Uft-off and consequently allows 
(^9) a lighting profile analysis along the entire ascent arc 

The angle ZALS, displayed in Fig. 16, is defined as the 
angle between the departure vector and the Sun-to-Earth 
(20) direction vector. It allows some judgment on available ascent 
lighting. 


Figure 16 displays a 3‘dimensional spatial view of the same 
typical launch geometry example shown previously in map 
format in Fig. 14. The difference in available range angles as 
well as orientation of the trajectory planes for the two daily 
launch opportunities clearly stands out. in addition, the figure 
illustrates the relationship between the “first” and “second 
daily” launch windows, defined in asymptote-relative time, 
^RLT' as contrasted with “morning” or “night” launches, 
defined in launch-site-local solar time. The latter is associated 


The length of the range angle required exhibits a complex 
behavior-the first launch window of the example in Figs. 14 
and 16 offers a longer range angle than the second, but the 
second launch window opens up with a range angle so short 
that direct ascent into orbit is barely possible. Further launch 
delay shortens the range even further, forcing the acceptance 
of a very long coast (one full additional revolution in parking 
orbit) before transplanetary departure injection. A detailed 
analysis of required arc lengths for the various sub-arcs of the 
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6^ BURNING ARC OF BOOSTER VEHICLES INTO PARKING ORBIT 

burning arc of final stage thrust 

$ RANGE ANGLE BETWEEN LAUNCH SITE AND DEPARTURE RADIAL ASYMPTOTE 

fio angle between perigee and departure radial asymptote 

ly TRUE ANOMALY OF INJECTION 

Fig. 17. Basic 9eomatry of the launch and ascent profile in the trajectory plane (after Ref. 9) 


departure trajectory h thus a trajectory design effort of para- r - perlapse radius, typicafly 6563 km, for a horizontal 
mount importance (Fig. 1 7). injection from a ] 85-km ( 100-nmi) parking orbit 


3. Range an^e arithinettc> For a viable ascent trajectory 
design, the range angle & must first of all accommodate the 
twin burn arcs and representing ascent into parking 
orbit and transplanetary injection bum into the departure 
hyperbola. In addition, it must also contain the angle from 
pcriapse to the direction, called “true anomaly of the 
asymptote direction’’ (Fig. 18): 


where: 


COSIL 


1 


3 p 


( 21 ) 


jj^ ~ GM, gravitational parameter of Earth (refer to the 
Table of Constants, Section V). 


The proper addition of these trajectory sub-arcs also requites 
adjustment for nonhorizontal injection (i.e„ for the flight 
path angle jf > 0), especially significant on direct ascent 
missions (no coast arc) and missions with relatively low 
dirust/weight ratio injection stages. The adjustment is accom- 
plished as follows: 


‘'i 


( 22 ) 


where: 
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Fig. 18. Angle from perigee to departure asymptote 


Vj = is the true anomaly of the injection point, usually 
is near 0 deg, and can be computed by iteration 
using: 


tan 7} 


1 + COS Vj 


(23) 


H 


injection fhght path angle above local horizontal, 
deg 

eccentricity of the departure hyperbola: 


= 1 + 




(24) 


To make Eq. (22) balance, the parking orbit coast arc, 
^COAST^ must pick up any slack remaining, as shown in 
Fig. 17. A negative ^coast that the -solution was 

too short-ranged. A direct ascent with positive injection true 
anomaly Vj (i.e., upward climbing flight path angle, jj, at 
injection) with attendant sizable gravity losses, may be accept- 
able, or even desirable (within limits) for such missions. 
Alternately, the other solution for 2^, exhibiting the longer 
range angle 6, and thus a longer parking orbit coast, 
should be implemented. An extra revolution in parking orbit 
may be a viable alternative. Other considerations, such as 
desire for a lightside launch and/or injection, tracking ship 
location and booster impact constraints, may all play a signif- 


icant role in the ascent orbit selection. A limit on maximum 
coast duration allowed (fuel boil-off, battery life, guidance 
gyro drift, etc.) may also influence the long/short parking orbit 
decision. In principle, any number of additional parking orbit 
revolutions is permissible. Shuttle launches of interplanetary 
missions (e.g., Galileo) are in fact required to use such addi- 
tional orbits for cargo bay door opening and payload deploy- 
ment sequences. In such cases, however, the precessional 
effects of Earth’s oblateness upon the parking orbit, primarily 
the regression of the orbital plane, must be considered. 

4. Parking orbit regression. The average regression of the 
nodes (i.e., the points of spacecraft passage through the 
equator plane) of a typical direct (prograde) circular parking 
orbit of 28.3-deg inclination with the Earth's equator, due to 
Earth's oblateness, amounts to about 0.46 deg of westward 
nodal motion per revolution and can be approximately com- 
puted from 


540° X X y X cos i 

n , deg/revolution (25) 

O 

where 

r = Earth equatorial surface radius, 6378 kr* 

= circular orbit radius, typically 6748 km for an 
orbital altitude of 370 km (200 nmi) 
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/j = 0.00108263 for Earth 

i - parking orbit inclination, deg. Can be computed for 
a given launch geometry from 

cos i = cos 0^ X sin 2^ (26) 

This correction, multiplied by the orbital stay time of N 
revolutions, must be considered in determining a biased launch 
time and, hence, the right ascension of the launch site at 
Uft-off: 

o, = a, + n X iV , deg (27) 

5. Dogleg ascent. Planar ascent has been considered exclu- 
sively, thus far. Reasons for performing a gradual powered 
plane change maneuver during ascent may be many. Inability 
to launch in a required azimuth direction because of launch 
site constraints is the prime reason for desiring a dogleg ascent 
profile. Other reasons may have to do with bum strategies 
or intercept of an existing orbiter by liie ascending spacecraft, 
especially if its inclination is less than the latitude of the 
launch site. Doglegs are usually accomplished by a sequence of 
out-of-plane yaw turns during first- and second-stage bum, 
optimized to minimize performance loss and commencing as 
soon as possible after the early, low-altitude, high aerodynamic 
pressure phase of flight is completed, or after the necessary 
lateral range angle offset has been achieved. 

By contrast, powered plane change maneuvers out of 
parking orbit or during transplanetary injection are much less 
efficient, as a much higher velocity vector must now be rotated 
through the same angle, but they may on occasion be opera- 
tionally preferable. 

As already discussed, a special geometric situation develops 
whenever the departure asymptote declination magnitude 
exceeds the latitude of the launch site, causing a “split azi- 
muth” daily launch window. Figure 13 shows the effects of 
asymptote declination and range safety constraints upon the 
launch problem. As the absolute value of declination increases, 
it eventually reaches the safety constraint on azimuth, prevent- 
ing any further pLnar launches. The situation occurs mostly 
early and/or late in the mission's departure launch period, and 
is frequeri^y associated with dual launches, when month-long 
departure periods are desired. 

The Shuttle-era S :e Transportation System (STS), includ- 
ing contemplated upper stages, is capable of executing dogleg 
operations as well. These would, however, effectively reduce 
the launch vehicle’s payload (or C 3 ) capability, as they did on 
expendable launch vehicles of the past. 


6 . Tradcuig and orientation. As the spacecraft moves away 
from the Earth along the asymptote, it is seen at a nearly 
constant declination- that of the departure asymptote, 6 ^ 
(DLA in the plotted data). The value of DLA greatly affects 
tracking coverage by stations located at various latitudes; 
highest daily spacecraft elevations, and thus best reception, are 
enjoyed by stations whose latitude is closest to DLA. Orbit 
determination, using radio doppler data, is adversely affected 
by DLA’s near zero degrees. 

The spacecraft orientation in the first few' weeks is often 
determined by a compromise between communication (antenna 
pointing) and solar heating constraints. The Sun-spacecraft- 
Earth (SPE) angle, defined as the angle between the outgoing 
F-infinity vector and the Sun-to-Earth direction, is very 
useful and is presented in the plots under the acronym ZALS. 
It was defined in the discussion of Fig. 16. This quantity has 
many uses: 

( 1 ) If the spacecraft is Sun-oriented, ^ALS equals the 
Earth cone angle (CA), depicted in Fig. 19 (the cone 



Fig. 19. Definition of cone and clock angle 
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angle of an object is a spacecraft-fixed coordinate, 
an angle between the vehicle's longitudinal -Z axis and 
the object direction). 

(2) The Sun-phase i ’igie, (phase angle is the Sun-object- 
spacecraft angle) ucscribes the state of the object’s 
disk lighting: a ^uUy lit disk is at zero phase. For the 
Earth (and Moon), several days after launch, the sun- 
phase ang.'e is: 

= 180-ZALS (28) 

(3) The contour labeled ZALS = 90'' separates two cate- 
gories of transfer trajectories- those early departures 
that first cut inside Earth’s orbit, thus starting out at 
negative heliocentric true anomalies for ZALS >90°, 
and those later ones that start ai positive true anoma- 
lies, heading out toward Jupiter and never experiencing 
the increased solar heating at distances of less than 
1 AU for ZALS <90°. 

7. Post4aunch spacecraft state. After the spacecraft has 
departed from the immediate vicinity of Earth (i.e., left the 
Earth’s sphere of influence of about 1-2 million km), it moves 
on a heliocentric conic, whose initial conditions may be 
approximated as 


^EARTH 

(29) 

\arth *** ^S/C ^ 

(30) 


where R and V of Earth are evaluated from an ephemeris at 
time of injection and Af represe* s time elapsed since then (in/* 
seconds). Th- vector in EME50 cartesian coordinates cwi 
be constructed using y/C^L = DLA = 6^, and RLA ='a^ 
in three components as follows: 

V = ( F X cos a X cos 5 , K X sin a X cos 6 , 

00^00 oo OQ ^ oe 00.00 

^'»Xsin5^) (31) 

8. Orbital launch problem. Orbital launch from the Shuttle, 
from other elements of the STS. or from aiiy temporary or 
permanent orbital space station complexes, introduces entirely 
new concepts into the Earth-departure problem. Some of 
the new constraints, already mentioned, limit our ability 
to launch a given interplanetary mission. The slowly regressing 
space station orbit (see Eq. 25) generally does not contain 
the K-infinity vector required at departure. Orbit lifetime 
or other considerations may dict'ite a space station’s orbital 
altitude that may be too bigli for an efficient injection burn. 
Innovative departure strategics are beginning to emerge, 
attempting to alleviate these problems - a recent Science 


Applications, Inc. (SAI) study (Ref. 10) points to some of 
the techniques available, such as passive wait for natural 
alignment of the continuously regressing space station orbit 
plane (driven by Earth’s oblateness) with the required 
K-infinity vector, or the utilization of 2- and 3- impulse man- 
euvers, seeking to perform spacecraft plane changes near the 
apogee of a phasing orbU where velocity is lowest and thus 
turning the orbit is easiest. These two approaches can be 
combined with each other, as well as with other suitable 
maneuvers, such as: 

(1) Deep space propulsive burns for orbit shaping and 
phasing, 

(2) Gravity assist flyby, including ;cturn AVF.GA, 

(3) Aerodynamic turns at grazi »^^ees or at imer- 

mediate planetary swingbys, ar. 

(4) Multiple revolution injection burns, requiring several 
low, grazing passes, combined with apogee plane change 
maneuvers, etc. 

All of these devices can be optimized to permit satisfactory 
orbital launches, as well as to achieve the most desirable condi- 
tions at the final arrival body. In general, space launch advan- 
tages, such as on-orbit assembly and checkout of payloads and 
clustered multiple propulsion stages, or orbital construction of 
bulky and fragile subsystems (solar panr s, sails, antennas, 
radiators, booms, etc.) will, it is hoped, greatly outweigh the 
significant deep-space mission penalties incurred because of 
the space station’s inherent orbital orientation incompatibility 
with departure requirements. 

D. Planetary Arrival Synthesis 

The planetary arrival trajectory design problem involves 
satisfying the project’s engineering and science objectives at 
the target body by shaping the arrival trajectory in a suitable 
Manner. As these objectives may be quite diverse, only four 
illustrative scenarios shall be discussed in this section -flyby, 
orbiter, atmospheric probe, and. to a small extent, lander 
missions. 

1. Flyby trajectory design. In this mission mode, the arrival 
trajectory is not modified in any deterministic way at the 
planel-the original aim point and arrival time are chosen to 
satisfy the largest number of potential objectives, long before- 
hand. 

This process involves the choice of arrival date to ensure 
desirable characteristics, such as the values of the variables 
VHP, DAP, ZAPS, etc., presented in plotted form in the data 
section of this volume. 
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Anothe. pair of significant variables on which to base 
arrival date selection are ZAPS and ZAPE-the angles between 
K-infinity and the planet-to-Sun and -Earth vectors, respec- 
tively. These two angles represent the cone angle (CA) of the 
planet during the far-encounter phase for a Sun- or an Earth- 
oriented spacecraft, in that order. ZAPS also determines the 
phase angle, of the planet’s solar illumination, as seen by 
the spacecraft on its far-encounter approach leg to the planet: 

= 180 -ZAPS (35) 

Both the cone angle and the phase angle have already been 
defined and discussed in the Earth departure section above. 

The flyby itself is specified by the aim point chosen upon 
ihe arrival planet target plane. This plane, often referred to 
as the ^-plane, is lighly useful aim point design tool. It is a 
plane passed through the center of a celestial body normal to 
. the relative spacecraft incoming veIoci^> vector at infinity . 
The incoming asymptote, i.e., the straight-line, zero-gravity 
extension of the vector, penetrates the ^-plane at the aim 
point. This point, defined by the target vector B in the ^-piane, 
is often described by its two components B • T and B • R, 
where the axes T and R form an orthogonal set with V^. The 
T-axis IS chosen to be parallel to a fundamental plane, usually 
the ecliptic (Fig. 21 ) or aiternatively . the planet’s equator. The 
magnitude of B equals the semi-minor axis of the flyby hyper- 
bola, b, and can be related to the closest approach distance, 
also referred to as the penapse radius, r^, by 

or 

''"n* (inf 

The direction angle ^ of the ^-vector, B, measured in the 
target plane clockwise from the T-axis to the ^-vector position 
can easily be related to the inclination, /, of the flyby trajec- 
tory, provided that both 6^ (DAP) and the T-axis, from which 
$ is measured clockwise, are defined with respect to the same 
fundamental plane to which the inclination is desired. For a 
system based on the planet equator (Fig. 22): 

cos = cos X cos (38) 

which assumes that is computed with the T-axis parallel 
to the planet equator (i.e., Tpp^ = X ot the 

ecliptic, as is frequently assumed (T£*^2: ” ^ 

Care must be taken to use the 8 angle as defined and intended. 



B MISS PARAMETER, BIS 
(TARGET VECTOR) 

0 AIM POINT ORIENTATION 
S PARALLEL TO INCOMING ASYMPTOTE, 
f PARALLEL TO ECLIPTIC PLANE 
AND i TO? 

R =SXT 

Fig. 21. Definition of target or arrival ^-plane coordinates 

The two systems of ^-plane T-axis definition can be recon- 
ciled by a planar rotation. -A0, between the ecliptic T- and 
R-axes and the T- and ft-axis orientations of the equator 
based sy stem 

sin (g^ - 

cos 6^ X tan ^EP - sin 6^ X cos (a„ - 

(39) 

where 

a^rp and ^^^p are right a^^cension and declination of the 
ecliptic pole in planet equatorial coordinates. For Venus 
using constants in Section V: 

= 155.4411, = 88.7250, deg 

0 ^ and 6^ are RAP and DAP, the directions of incoming 
V ’ m pia. equatorial coordinates. 
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ECLIPTIC POLE (EP) 



^PEQ “ ^ ^PEQ 

Aor = 90 - [(a^p - 90) - (a„ - 180)] =a^- a^p 

A0= ARCTAN [SIN Aa/(COS 5^ x TAN b^p - SIN 6^ x COS Aa)l 

Fig. 22. Two Taxis definitions in the arrival ^>plane 


The correction A0 is applied to a 0 angle computed in the 
ecliptic system as foMows(Fig. 22): 


^PEQ 




(40) 


The ecliptic T, axes, however, have to be rotated by 
(clockwise direction is positive in the ^-plane) to obtain 
planet equatorial T, coordinate axes The ^-magnitude 

of an aim point in either system is the same. 


periapse, at the entry point of a probe, or generally at any 
position r (subscript S - Sun, could be replaced by E- Earth, 
if desired), (see Fig. 24): 

cos = -cos X cos ZAP^ - sin X sin ZAP^ 

Xco%{ET^P-6^^^) (41) 

where 


The projections of the Sun-to-planet and Earth-to-planet 
vectors into the £^p^ane represent aim point loci of diametric 
Sun and Earth occultations, respectively, as defined in Fig. 23. 
The ^-plane 0-angles (with respect to Tecz, axis) of these vari- 
ables are presented and labeled ETSP and ETEP, respectively, 
in the plotted mission data. In addition to helping desi, n or 
else avoid diametric occultations, these quantities allow 
computation of phase angles, of the planet at the spacecraft 


^SjC “ point angle in the ^-plane, must be with 

respect to the same T, as ETSP. 

= the arrival range angle from infinity (a position far 
out on the incoming asymptote) to the point of 
interest r (Fig. 25). 

The computation of the arrival range angle, to the 
position of the desired event depends on its type, as follows: 
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where is the true anomaly at entry, should always 
be negative, and can be computed if entry radius and 
altitude, and 7^^, the entry angle, are 

known: 



whereas the fictitious periapse radius for the entry, 
to satisfy 7f at is equal to 



(47) 


The B value corresponding to this entry point can be com- 
puted from Eq. (36), while the 8 angle in the ^plane would 
depend on the desired entry latitude inclination, Eq. (38), or 
phase angle, Eq. (41). 

The general flyby problem poses the least stringent con- 
straints on a planetary encounter mission, thus allowing 
optimization choices from a large list of secondary parameters, 
such as satellite viewing and occultation, planetary fields and 
particle in situ measurements, special phase-angle effects, etc. 

A review of the plotted handbook variables, required in the 
phase-angle equation (Eq. 41 ), shows that the greatest magni- 
tude variations are experienced oy the ZAPS angle, which is 
strongly flight-time dependent: the longer the trip, the smaller 
ZAPS. For low equatorial inclination, direct flyby orbits, this 
implies a steady move of the periapse towards the lit side and, 
eventually, to nearly subsolar periapses for long missions. This 
also implies that on such flights the approach legs of the tra- 
jectory are facing the morning terminator or even the dark 
side, as trip time becomes longer, exhibiting large phase angles 
(recall that phase is the supplement of the ZAPS angle on 
the approach leg). This important variation is caused by a grad- 
ual shift of the incoming approach direction, as flight time in- 
creases, from the subsolar part of the target planet's leading 
hemisphere (in the sense of its orbital motion) to itsantisolar 
part- 

The arrival time choice on a very fine scale may greatly 
depend on the desire to observe specific atmospheric/surface 
features or to achieve close encounters with specific satellite^ 
of the arrival planet. Passages through special satellite event 
zones, e.g., flux tubes, wakes, geocentric and/or heliocentric 


occultations, require close control of arrival time. The number 
of satellites passed at various distances also depends on the 
time of planet CjA. These fine adjustments do, however, 
demand arrival time accuracies substantially in exc'ssof those 
provided by the computational algorithm used in if ef’“ r» 
which generated the subject data (accuracies of 1-5 min for 
events or 1-2 h for encounters would be required vs uncertain- 
ties of up to 1.5 days actually obtained with the rectilinear 
impact pseudo-state theorem). Numerically searched-in inte- 
grated trajectories, based on the information presented as a 
first guess input, are mandatory for such precision trajectory 
work. 

Prehni inary design considerations for penetrating, grazing, 
or avoiding a host of p\:net v e».tered fields \ ,*:;icle struc- 
tures. such as magnetic fields, radiation belts, plasma ton, ring 
and debris structures, Owwultations by Sun, Earth, stars, or 
satellites, etc., can ail be presented on specialized plots, e.g.. 
the i^plane. and do affect the choice of suitable aim point 
and arrival lime. All of these studies require the propagation 
of a number of flyby trajectories. Adequate initial conditions 
for such efforts can be found in the handbook as: VHP ( 
and DAP (6^) already defined, as well as RAP (ot^), the 
planet equatorial right ascension of the incoming asymptote 
(i.e.. its east longitude from the ascending node of the planet's 
mean orbital plane on its mean equator, both of date). The 
designer's choice of the aim point vector, either as £ and d, or 
as cartesian B • T and B * R, completes i* e input set. Suitable 
programs generally exist to process this information. 

2. Capture orbit design. The capture problem usually 
involves the task of determining what kind of spacecraft orbit 
is most desired and the interconnected problem of how and at 
what cost such an orbit may be achieved. A scale of varying 
complexity may be associated with the effort envisioned-an 
elliptical long period orbit with no specific orientation at the 
trivial end of the scale, through orbits of controlled or opti- 
mized lines of apsides (i.e., periapse location), nodes, inclina- 
tion, or a safe perturbed orbital altitude. Satellite G/A-aided 
capture, followed by a satellite tour, involving multiple satel- 
lite G/A encounters on a number of revolutions, each designed 
to achieve specific goals, probably rates as the most complex 
capture orbit class. Some orbits are energetically very difficult 
*0 achieve, such as close circular orbits, but all require signifi- 
cant expenditures of fuel. As maneuvers form the background 
to this subject a number of useful orbit design concepts shall 
be presented to enable even an unprepared user to experiment 
with the data presented. 

The simplest and most efficient mode of orbit injection is 
a coplanar burn at a common periapse of the arrival hyperbola 
and the resulting capture orbit (Fig. 26). The maneuver A V 
required is: 
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I — T ,km/s (48) 


A plot of orbit insertion AV required as a function of and 
P (using Eq. 50) is presented in Fig. 27. The apoapse radius of 
such an orbit of given period would be 


The orbital period for such an orbit, requiring knowledge 
of periapse and apoapse radii, and , is 



If on the other hand, a known orbit period P (in seconds) 
is desired, the expression for AFis 



2^XP^ 

; tp .km (51) 

An evaluation of Eq. (50) (and Fig. 27) shows that lowest 
orbit insertion AF is obtained for the lowest value of tlie 
longest period F, and the lowest of arrival. 

Of some interest is injection into circular capture orbits, a 
special case of the coapsidal insertion problem. It can be 
shown (Ref. 9) that an optimal AF exists for insertion into 
capture orbits of constant eccentricity, including e = 0, i.e., 
circular orbits, which would require a specific radius: 




Fig. 26. Coapsidal and corangential capture orbit insertion geometries 
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Fig. 27. Coapsidal capture orbit insertion maneuver AV 
requirements for Venus (using Eq. 50) 


’CO 



(52) 


while the corresponding optimal value for A V would be 




CO 



(53) 


Frequently the orbital radius obtained by use of Eq. (52) is 
incompatible with practical injection aspects or with arrival 
planet science and engineering objectives. 

A more general coplanar mode of capture orbit inseition, 
requiring only tangentiality of the two trajectories at f»n 
arbitrary maneuver point of radius r common to both orbits. 
Fig. 26. requires a propulsive effort of 


AV 



2 fji {r. ^ “ '*) 

'■(»•, + >• ) 

/I p' 


(54) 


It can be clearly seen that by performing the burn at periapse 
the substitution r = brings us back to Eq. (48). 

The cotangential maneuver mode provides nonoptimal con- 
trol over the orientation of th * major axis of the capture orbit. 
If it is desired to rotate this line of apsides clockwise by 


one can solve for the hyperbolic periapse and the 
bum radius r using selected values of true anomaly at hyper- 
bolic bum point and its capture orbit equivalent 

"f “ "h + (55) 

Utilizing the following three equations (where E and H stand 
for elliptic and hyperbolic, respectively): 



The procedure of obtaining a solution to these equations is 
iterative. For a set of given values for r,, r , and an assumed 

A p 

Aco, a set of and the hyperbolic and elliptical burn 
point true anomalies which would satisfy Eqs. (55-58) can 
be found. This in turn leads to r, the maneuver point radial 
distance, and hence, AV (Eq. 54). A plot of AK cost for a set 
of consecutive Acu^ choices will provide the lowest AV value 
for this maneuver mode. 

For an optimal insertion into an orbit of an arbitrary major 
axis orientation one must turn to the more general, still 
planar, but intersecting (i.e., nontangential burn point) man- 
euver (see Fig. 28). It provides sufficient flexibility to allow 
numerical optimization of AV with respect to apsidal rotation, 

A more appropriate way to define apsidal orientation is to 
measure the post-maneuver capture orbit periapse position 
angle with respect to a fixed direction, e.g., a far encounter 
point on the incoming asymptote, thus defining a cap- 
ture orbit periapse range angle: 
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1ST BURN POINT 
OPTION 


PROVIDED IS AVOIDED. 

OPTIMAL BURN AV 
REQUIRES 


Fig. 28. Coapsidal and intersecting capture orbit insertion geometries 


A(jO = - 

OO p oo 



Taken from Ref. 9, the expression for the intersecting burn 
Al^is: 



where 

« y t^p ('• - ''cv. ) + ‘ « <'■' - ^Ia )i (•- - '•/-) K - '■> j 

(60) 

r = planet centered radius at burn, 
km 

and = apoapse and periapse radii of capture ellipse 
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tcA = closest approach radius of flyby hyperbola 

6 - flag: 6 = +1 if injection occurs on same leg 

(inbound or outbound) of both hyperbola 
and capture ellipse, 6 = -1 if not. 

It should be pointed out that Eq. (57) and (58) still apply 
in the intersecting insertion case, while Eq. (56) does not 
(as it assumes orbit tangency at burn point). 

The evaluation of intersecting orbit insertion is more 
straightforward than it was for the cotangential case. By 
assuming and orbit size (e.g.. = 4 km's, = 1.04% 

RV,P= 24 hours, a typical capture orbit) and stepping thiough 
a set of values for the capture orbit burn point true anomaly, 
one obtains, using Eqs. (57-60), a family of curves, one for 
each value of hyperbolic closest approach distance. 

As shown in Fig. 29, the envelope of these curves provides the 
optimal insertion bum AK for any value of apsidal rotation 
Aco^ desired. The plot also shows clearly that cotangential and 
apsidal insertion burns are energetically inferior to burns on 
the envelope locus. For the same assumed capture orbit, a 
family of optimal insertion envelopes, for a range of values 
of arrival is presented in Fig. 30. 

The location of periapse and node with respect to the sub- 
solar point is of extreme importance to man> mission objec- 
tives. It can be controlled by choice of departure and arri\al 
dates, by Ar expenditure at capture orbit insertion, by an 
aerodynamic maneuver during aerobraking. b\ depending on 
the planet's motion around the Sun to move the subsolar point 
in a manner optimizing orbital science, or by using natural 
perturbations and making a judicious choice of orbit size, 
inclination, i. and initial argument of periapsis. such as to 
cause regression of the node, fi, and the advance of periapsis. 
ci, both due to oblateness and solar attraction, to move the 
orbit m a desired manner or at a specific rate. For an elliptical 
capture orbit (Fig. 31), the rates due to oblateness (with /, 
n, and CO referenced to the equator plane) are. 

• - 3 '^1 180 

= — X ^ — cos i X , deg/s (61) 

00=4^ (2-(5/2)sin^ /)X — ,deg/s 

p2 

(62) 

where 

n = mean orbital motion, rad/s 

a = (r^ + rp)/2. semi-major axis of capture orbit, km 
p = 2r^ rp)' semi-Iatus rectum, km 


= Equatorial surface radius of \ enus, km 

J 2 - Oblateness second harmonic coefficient of Venus (for 
value see Sv :tion V on constants). 

It should be noted that = 0 occuis for / = 90 deg, while 
CO = 0 is found for / - 63.435 deg. 

For the grazing 24-hour Venus orbiter example of Figs. 29 
and 30, maximum values for Q and co (-.0006 and .001 1 degy 
day, respectively) are negligible, due to the very small value of 
J 2 at Venus. Even for a grazing circular orbit (h = 300 km) the 
two rates would amount to only -.034 and .068 deg/day. For 
the same reason, the effects of the higher harmonics. and 
would produce change rales of similar order with those due 
to the small Because of the relative proximity of 

Venus to the Sun. solar third-body effects upon the 24-hour 
grazing orbit should be considered, amounting up to 17 = 
-.033; cj would be of the same order. Solar perturbations of 
periapse radius, are more critical and could amount up to 
about ±4.5 km/day for the same sample orbit. From Refs. 1 1- 
13: 



where the spacecraft ciemerts / and cc are now referenced to 
the plane of Venus' orbit (about the Sun), while and 
describe Venus' mean motion and eccentricity, = 27 t, 
(224.7007 X 86400. ) and 6y ^ 0. 

3. Entry probe and lander trajectory design. Entry tra- 
jectory design is on one hand concerned with maintenance of 
acceptable probe entry angles and low relative velocity with 
respect to the rotating atmosphere. On the other hand, the 
geometric relationship of entry point, subsolar point and 
Earth (or relay spacecraft) is of paramount importance. 

Lighting during entry and descent is often considered the 
primary problem to be resolved. As detailed in the flyby and 
orbital sections above, the choice of trip time affects the value 
of the ZAPS angle which in turn moves the entry point for 
longer missions closer to the subsolar point and even beyond, 
towards the morning terminator. 

Landers or balloons, regardless of deceleration mode, 
prefer the morning terminator entry point which provides a 
better chance for vapor-humidity experiments, and allows a 
longer daylight interval for operations following arrival. 

The radio-link problem, allowing data flow directly to 
Earth, or via another spacecraft in a relay role, is very complex. 
It could require studies of the Earth phase angle at the entry 
locations or alternately, it could require detailed parametric 
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Fig. 31. General satellite orbit parameters (from Ref. 9) 


studies involving relative motions of probe and relay spacecraft 
throughout probe entry and its following slow descent. 

BaUoon missions cuuld also involve consideration oi a 
variety of wind drift models, and thus, are even more complex 
as far as the communications problem with the Earth or the 
spacecraft is concerned. 

E. Launch Strategy Construction 

The constraints and desires, briefly discussed above, may be 
displayed on the mission space launch/arrival day plot as 
being limited by the contour boundaries of C^L, the dates, 
DLA, VHP, ZAP, etc., thus displaying the allowable launch 
space. 

Within this launch space a preferrt^d day-by-day launch 
strategy must be specified, in accordance with prevailing 
objectives. The simplest launch strategy, often used to maintain 
a constant arrival date at the target planet, results in daily 
launch points on a horizontal line from leftmost to rightmost 
maximum allowable C^L boundary for that arrival date. 
Such a strategy makes use of the fact that most arrival charac- 
teristics may stay nearly constant across the launch space. 
Lighting and satellite positions in this case are fixed, thus 
allowing a similar encounter, satellite G/A, or satellite tour. 


A different choice of strategy could be to follow a contour 
line of some characteristic, such as DLA or ZAP. One could 
also follow the minimum value locus of a parameter, e.g., 
C 3 Z (i.e., the boundary between Class 1 and 2 will tin Type 1 
or II) for each launch date, throughout the launch space. 

Fundamentally different is a launch .strategy for a dual or 
multiple spacecraft mission, involving more than one launch, 
either of which may possibly pursue divergent objectives. As 
an example. Fig. 32 shows the Voyagers 1 and 2 launch strat- 
egy, plotted on an Earth departure vs Saturn arrival date plot. 
A 14-day pad turnaround separation between launches was 
to be maintained, a 10 -day opportunity was to be available 
for each launch, and the two spacecraft had substantially 
different objectives at Jupiter and Satum-one was to be 
lo-intensive and a close Jupiter flyby, to be followed by a 
close Titan encounter at Saturn, and the other was Ganymede- 
and/or Callisto-intensive, a distant Jupiter flyby, a' a safety 
precaution against Jovia n radiation damage, aimed to continue 
past Saturn to Urami' ' * Neptune. Here, even t* ^ spacecraft 
departure order wa. .versed by the strategy within the 
launch space. 

Launch strategies for orbital departures from a space station 
in a specific orbit promise to introduce new dimensions into 
mission p nning and design. New concepts are beginning to 
emerge on this subject e.g.. Refs. 10 and 14. 


OCALLISTO INTENSIVE # GANYMEDE AND 10 iNTE'iSIVt 



Fig. 32. Voyager (MJS77) trajectory space and launch strategy 
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IV. Description of Trajectory Characteristics 
Data 

A. General 

The data represent tiajectory performance information 
plotted in the departure date vs arrival date space, thus defln> 
ing all possible direct ballistic transfer trajectories between the 
two bodies within the time span considered for each oppor- 
tunity. Twelve individual parameters are contour-plotted. The 
first, C 3 /., is plotted bold on a Time of Flight (TFL) back- 
ground; the remaining ten variables are plotted with bold con- 
touring on a faint C^L background. Eleven plots are presented 
for each of eight mission opportunities between 1991 and 
2005. 

The individual pl^ts are labeled in the upper outer corner 
by bold logos displaying an acronym of the variable plotted, 
the mission's departure year, ar : a symbol of the target planet. 
T ese permit a quick and fail-safe location of desired informa- 
tion. 

B. Definition of Departure Variables 

C^L\ Earth departure energy (km^/s^); same as the 
square of departure hyperbolic excess velocity 
Vl = C^L=V] -2 %//?,. where 

= conic injection velocity (km/s). 

injection radius (km), sum of 
surface radius planet injection 
altitude hj, where 

Earth's surface radius. (For value, see 
Section V on constants.) 

= gravitational constant times mass of the 
launch body (for values, refer to Section V 
on constants). 

C^L must be equal to or exceeded by the launch 
vehicle capabilities. 

DLA: , geocentric declination (vs mean Earth equator 

of 1950.0) of the departure vector. May im- 
pose laundt constraints (deg). 

RLA: , geocentric right ascension (vs mean Earth 

equator and equinox of 1950.0) of the departure 
vector. Can be used with C^L and DLA to 
compute a heliocentric initial state for trajectory 
analysis (deg). 

ZALS: Angle between departure vector and Sun-Earth 
vector. Equivalent to Earth-p >be-Sun angle several 
days out (deg). 

C. Definition of Arrival Variables 

VHP. planetocentric arrival hyperbolic excess 

velocity or F-inflnity (km/s), the magnitude of the 


vector obtained by vectorial subtraction of the 
heliocentric planetary orbital velocity from the 
spacecraft arrr'al heliocentric velocity. It repre- 
sents planet-relative velocity at great distance from 
target planet, at beginning of far encounter. Can 
be used to compute spacecraft velocity at any 
point r of Hvby, including CjA (periapse) d*s- 
tance : 



where 

MP(VENUS) “ gravitational parameter GM of the arrival 
planet. (For the value, refer to Section V 
on constants.) 

DAP: 5oo^ , planetocentric declination (vs mean planet 
equator of date) of arrival vector. Defines 
lowest possible flyby/orbiter equatorial inclination 
(deg) 

RAP , planetocentric right ascension (vs mean 
planet equator and equinox of date, i.e., RAP is 
measured in the planet equator plane from ascen^ 
ing node of the planet's mean orbit plane on the 
planetary equator, both of date). Can be used 
together with VHP and DAP to compute an 
initial flyby trajectory state, but requires ^-plane 
aim point information, e.g., B and B (deg). 

ZAPS: Angle between arrival vector and the arrival 
planet-to-Sun vector. Equivalent to planet-probe- 
Sun angle at far encounter; for subsolar impact 
would be equal to 180 deg. Can be used with 
ETSP, VHP, DAP, and B to determine sola. ,^hase 
angle at periapse, entry, etc. (deg). 

ZAPE: Angle between arrival vector and the planet-to- 
Earth vector. Equivalent to ^^anet-probe-Earth 
angle at far encounter (deg). 

ETSP: Angle in arrival 5 -plane, measured from T-axis*, 

clockwise to projection of Sun-to-planet vector. 
Equivalent to solar occultation region centerline 
direction in 5-plane (deg). 

ETEP: Angle in arrival 5 -plane, measured <Vom T-axis*, 
clockwise, to projection of Earth-to-planet vector. 
Equivalent to Earth occultation region centerline 
direction in 5-plane (deg). 


♦ETSP and ETEP plots arc based on T-axis defined as being parallel to 
ecliptic plane (see text for explanation). 
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-1.3 X 10“^ 


V. Table of Constants ^ 

Constants used to generate the information presented aie j - -2 4 X 10 ^ 

summarized in this sect ion . v en us 


A. Sun 

GM = 132,712,439,935. km^/s^ 
= 696,000. km 


Solid surface radius: 

/?venus 6051.813 km 

SURFACE 

Direction of the Venus planetary equatorial lAU north 
pole (in Earth Mean Equator of 1950.0 coordinates). 


B. Earth/Moon System 

GJUsystem ^ 403,503.253 km^/s^ 
GM^arth " 398.600.448 073 km^/s^ 


0 ^ = 272.80 deg. 8^= 67.20 deg 
Siderial rotation period (retrograde) 


= 0.00108263 

INEARTH = km 

SURFACE 


C. Venus 

^^VENUS " 324,858.7656 km 


h 


VENUS 


4.5 X 10“^ 


P - 243.010 days 

D. Sources 

The constants represent the DE-118 planetary ephemeris 
(Ref. 15) and Maimer lO/Piofieer Venus Orbiter trajectory 
reconstruction da;a. Definition of th^ Earth’s equator 
(EME50.0) is consistent with Refs. 16-17, but would require 
minor adjustments for the nt.w equator and equinox, epoch 
of J2000.0 (Ref. 18). Definition of the Venus po^e is also in 
accordance with Ref. 18. 
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